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ABSTRACT 

Cosmic magnification is due to the weak gravitational lensing of sources in tlie dis- 
tant Universe by foreground large-scale structure leading to coherent changes in the 
observed number density of the background sources. Depending on the slope of the 
background source number counts, cosmic magnification causes a correlation between 
the background and foreground galaxies, which is unexpected in the absence of lensing 
if the two populations are spatially disjoint. Previous attempts using submillimetre 
(sub-mm) sources have been hampered by small number statistics. The large num- 
ber of sources detected in the Herschel Multi-tiered Extra-galactic Survey (HerMES) 
Lockman-SWIRE field enables us to carry out the first robust study of the cross- 
correlation between sub-mm sources and sources at lower redshifts. Using ancillary 
data we compile two low-redshift samples from SDSS and SWIRE with (z) ~ 0.2 
and 0.4, respectively, and cross-correlate with two sub-mm samples based on flux 
density and colour criteria, selecting galaxies preferentially at z ~ 2. We detect cross- 
correlation on angular scales between ~ 1 and 50 arcmin and find clear evidence that 
this is primarily due to cosmic magnification. A small, but non-negligible signal from 
intrinsic clustering is likely to be present due to the tails of the redshift distribution 
of the sub-mm sources overlapping with those of the foreground samples. 
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1 INTRODUCTION 

Large-scale structure at low redshifts systematically magni- 
fies sources at higher redshifts as a result of gravitational 
light deflection in the weak limit. On the one hand, fewer 
sources will be observed, because lensing stretches the solid 
angle and dilutes the surface density of sources. Conversely, 
the effective flux limit is lowered as a result of magniflcation, 
which leads to a deeper survey. Whether there is an increase 
or decrease in the observed number density of sources de- 
pends on the shape of the background source number counts 
- an effect known as the magnification bias (Bartelmann 
& Schneider 2001; hereafter BSOl). At submillimetre (sub- 
mm) wavelengths the magnification bias is expected to be 
large and positive, resulting in an increase in the observed 
number density of sources compared to the case without 
lensing (e.g. Blain & Longair 1993; Blain et al. 1996; Ne- 
grello et al. 2007; Lima et al 2010). 

Cosmic magnification also induces an apparent angu- 
lar cross-correlation between two source populations with 
disjoint spatial distributions. It can thus be measured 
by cross-correlating non-overlapping foreground and back- 
ground samples. When combined with number counts, such 
a cross-correlation can provide constraints on cosmological 
parameters (e.g. fim, cs) and galaxy bias, a key ingredient in 
galaxy formation and evolution models (Menard & Bartel- 
mann 2002). As the weak lensing-induced cross-correlation 
also probes the dark matter distribution, it provides an in- 
dependent cross-check of the cosmic shear measurements, 
which depend on the fundamental assumption that galaxy 
ellipticities are intrinsically uncorrclatod. Most previous in- 
vestigations, using foreground galaxies selected in the opti- 
cal or infrared together with background quasars, have pro- 
duced controversial or inconclusive results (e.g. Seldner & 
Peebles 1979; Bartelmann & Schneider 1994; Bartsch et al. 
1997). The best detection to date is presented in Scranton 
et al. (2005), where cosmic magnification is detected at an 
8a significance level using 13 million galaxies and ~200,000 
quasars from the Sloan Digital Sky Survey (SDSS). 

The amplitude of the weak lensing-induced cross- 
correlation is determined by several factors: the dark matter 
power spectrum and growth function, the shape of the back- 
ground source number counts and the bias of the foreground 
sources. At sub-mm wavelengths, the power-law slope of the 
cumulative number count is exceptionally steep, >2.5 for 
sources in the flux range 0.02 — 0.5 ,Jy at 250, 350 and 
500 iim (e.g. Patanchon et al. 2009; Oliver et al. 2010a; 
Glenn et al. 2010; Clements et al. 2010). In Scranton et 
al. (2005), the number count slope of the quasar sample is 
considerably flatter (~ 2 for the brightest ones). In addi- 
tion, sub-mm sources detected in deep surveys mainly re- 
side in the high-redshift Universe with a median redshift of 
« ~ 2 (Chapman et al. 2003, 2005; Pope et al. 2006; Aretx- 
aga et al. 2007; Amblard et al. 2010). The steep number 



counts, together with the large redshift range, make sub-mm 
sources an ideal background sample. So far there have been 
two attempts at measuring the weak lensing-induced cross- 
correlation between foreground optical galaxies and back- 
ground sub-mm sources, but with conflicting results. Al- 
maini et al. (2005) measured the cross-correlation between 
39 SCUBA sources and optical sources at lower redshifts 
(z) ~ 0.5. They claimed evidence for a significant signal 
which might be caused by lensing. Conversely, Blake et al. 
(2006) did not find evidence for cross-correlation due to cos- 
mic magniflcation using a similar number of sources. 

The Herschel Multi-tiered Extra-galactic Survey (Her- 
MES, Oliver et al. 2010b) is the largest project being under- 
taken by Herschel (Pilbratt et al. 2010). In this paper, we 
calculate the angular cross-correlation between foreground 
galaxies selected from SDSS or the Spitzer Wide-area In- 
frared Extragalactic (SWIRE; Lonsdale 2003, 2004) survey 
and background sub-mm sources detected by the Spectral 
and Photometric Imaging Receiver (SPIRE; Griffin et al. 
2010) instrument on Herschel. This paper is organised as 
follows: In Section 2, we give a brief introduction to mag- 
nification bias and the angular cross-correlation function. 
In Section 3, we describe the various data-sets used as fore- 
ground and background samples. Measurements of the cross- 
correlation between foreground and background samples are 
presented in Section 4. Finally, discussions and conclusions 
arc given in Section 5. Throughout the paper, we use a spa- 
tially flat ACDM cosmology with fim = 0.3 and Ho — 70 
km s~^ Mpc~^. Magnitudes are in the AB system. 



2 MODELLING THE CROSS-CORRELATION 
FUNCTION 

In this section we briefly describe the magnification bias 
and how it manifests itself in the number density and cross- 
correlation between two spatially separated populations. We 
refer the reader to Moessner & Jain (1998), BSOl, Cooray & 
Sheth (2002) and references therein for a complete introduc- 
tion. Suppose a background population has an intrinsic (i.e. 
unlensed) number density nu{S,z), where S is fiux density 
and z is redshift. As a result of lensing, the sky solid angle 
is stretched locally by a factor of /x((/), z) {4> denotes angular 
position on the sky) and S is magnified by the same factor 
because surface brightness is preserved. The two contrasting 
effects modify the observed (lensed) number density in the 
following way 

ni(g,z) = ""^^/^('^'")'"^ (1) 

ti{4>,z) 

When the lens plane is at a much lower redshift than the 
source plane, the redshift-dependent magnification can be 
substituted by the magnification ^ of a source at infinity. 
Assuming the cumulative number count distribution of the 
background population can be described by a power-law 
Nu{S) oc , we should expect a factor of 



Detection of cosmic magnification of suh-mm galaxies 3 



(2) 



Ni{S) 

change in the observed number count. Strictly speaking, the 
number count slope /3 = /3(S) is a function of flux density. 
In this paper, we make the simplifying assumption that /3 is 
a constant over the flux range we probe. Using the number 
counts of resolved sources presented in Oliver et al. (2010a), 
we find that in the flux range 0.03 — 0.5 Jy, /3 = 2.53 ± 
0.16, 2.99 ± 0.51 and 2.66 ± 0.24 at 250, 350 and 500 ^m 
respectively. 

The angular cross-correlation function between popula- 
tion 1 at lower redshifts and population 2 at higher redshifts 
is defined as 



(3) 



where Srii = Ui (0) /fii — 1 is the number density fluctuation 
and fii is the average number density of the ith sample. We 
can decompose Srii into two parts. 



5n.(0) = (5nK0) + 5nf(0). 



(4) 



The flrst term Sn^ is due to intrinsic clustering of galaxies 
and is a projection of density fluctuations along the line-of- 
sight. 



(5) 



= / dxW,{x)5{r{x)^,a), 
Jo 

where xh is the comoving radial distance to the horizon, 
r(x) is the comoving angular diameter distance, W{x) is the 
normalised radial distribution of the sources, a is the scale 
factor, S{r{x)4'jO-) is the dark matter density perturbations 
and bi is the bias factor assumed to be scale- and time- 
independent. The second term in equation (4) Sn^ is due to 
magniflcation bias, 

Ni - Nu 



1 = 2(/3 - 1)k. 



(6) 



In the last step, we have used the weak lensing limit, /i = 
1 -I- 2k. The convergence k is a weighted projection of the 
density field along the line of sight (BSOl), 



^ Jo 

where g{x) is the radial weighting function defined as 
g(x) = r(x) ^^^^^W(x')dx'. 



(7) 



(8) 



The angular cross-correlation between the two populations 
is then 

+ {Sn>^{4>)Sn^{4:')) + {Sn'i{4>)SnU$')). (9) 

The first term l^5n\[cj>)5n2{(j>')) is due to the intrinsic clus- 
tering of the two populations tracing the same large-scale 
structure, 



2^ 



P{X,k)Jo{kre)dk, (10) 



lOcc(e) = bib2 / WiWidx 

Jo Jo 

where -P(Xi k) is the dark matter power spectrum and 
Jo (a;) = sin(a;)/a; is the zeroth-order Bessel function. Note 
that Wcc vanishes if the two populations have disjoint spatial 
distribution. The second term {Sni{(l})5n2 {<j>')} is caused by 




Figure 1. Redshift distribution of the foreground and back- 
ground populations normalised so that the peak of each N(z) 
is equal to unity. For the foreground sample Fl, the N(z) is de- 
rived from spectroscopic redshifts. For the foreground sample F2, 
we have either spectroscopic redshifts or good quality photomet- 
ric redshifts. The N{z) for the two background samples, Bl and 
B2, are derived from sub-mm colours using modified black-body 
spectra. 
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Wi^dx 
a 



— P{x,k)Jo{kre)dk. 
2n 



(11) 



The third term {5n'^{<^)5n2(4>')) is due to weak lensing by 
large-scale structure in front of both the foreground and 
background sources. The last term {5n'^{ij))5n2{4>')) repre- 
sents large-scale structure traced by the background sources 
lensing the foreground sources which is only present if the 
two samples have overlapping redshift distributions. The 
last two terms are negligible. To derive the expected cross- 
correlations (wcc and uifb) between our foreground and back- 
ground samples in Section 4, we use the CAMB software 
package (Lewis, Challinor & Lasenby 2000), which is based 
on CMBFAST (Seljak & Zaldarriaga 1996), to generate the 
non-linear matter power spectrum using the fitting formulae 
of Smith et al. (2003). 



3 DATA-SETS 

For the first foreground sample, referred to as Fl, we se- 
lect 7,761 sources with r < 19.4 from the SDSS DR7 in 
Lockman-SWIRE observed by Herschel-'SPIR'E. The star- 
galaxy separation is done in the same way as in Stoughton 
et al. (2002). The redshift distribution N{z) of the sample 
Fl is derived from spectroscopic redshifts obtained in the 
Galaxy and Mass Assembly (GAMA) survey (Baldry et al. 
2010). The median redshift of Fl is ~ 0.2. The second fore- 
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Table 1. Summary of foreground and background samples. The 
columns are the sample name, the number of sources, the median 
redshift and the selection criterion. For the two background sam- 
ples, Bl and B2, we list the number of sources in the overlapping 
region with Fl and F2, respectively. 



Sample 






Selection criterion 


Fl 


7,761 


- 0.2 


r < 19.4 


F2 


13,888 


~ 0.4 


53.6 > 100 ^liy 


Bl 


2,477 / 1,886 


- 2.0 


S350 > 35 mjy 


B2 


2,398 / 1,848 


~ 2.5 


S350/5'250 > 0.85 



ground sample F2 is selected from sources detected by the 
Spitzer Infrared Array Camera (IRAC; Fazio et al. 2004) in 
the SWIRE survey. Full details of the data processing and 
catalogues can be found in Surace et al. (2005). We select 
13,888 sources with S^.e ^ 100 fiJy in the region overlapping 
with the Herschel-SPTRE observation in Lockman-SWIRE. 
The star-galaxy separation is performed in the same way 
as in Waddington et al. (2007). About 17 percent of the 
sources in F2 have spectroscopic redshifts and ~80 percent 
have good quality photometric redshifts with redshift un- 
certainty < 0.1 (Oyaizu et al. 2008; Rowan- Robinson et al. 
2008). The median redshift of F2 is z ~ 0.4. 

To construct the background samples in the 13.2 deg^ 
Lockman-SWIRE field, we use the single-band SPIRE cat- 
alogues generated by the SUSSEXtractor source extractor 
in HIPE (Smith et al. 2010). The cross-match between the 
250 and 350 fim catalogue is done by selecting the bright- 
est 250 fim source within 12.6 arcsec of a 350 /im source 
(FWHM=25.2 arcsec at 350 fim). The flux density at which 
the integral source counts reach 1 source per 40 beams is 
18.7 and 18.4 mJy at 250 and 350 fim, respectively (Oliver 
et al. 2010a). The first background sample, Bl, comprises 
sources brighter than 35 mJy at 350 fim. In total, there are 
2,477 / 1,886 such sources in the region that overlaps with 
Fl / F2. The second background sample B2 includes sources 
with S-i5o/S250 > 0.85. In total, there are 2,398 / 1,848 such 
sources in the overlapping region with Fl / F2. About 50% 
of the sources in Bl are found in B2 as well. Because most of 
the background sources do not have spectroscopic redshifts, 
we make use of the sub-mm colours and modified black- 
body templates to generate qualitative redshift distributions 
which are consistent with typical model predictions (e.g. Le 
Borgne et al. 2009; Valiante et al. 2009). The majority of 
the sources with 5*350 ^ 35 mJy lie at 1.5 < z < 3 and peak 
at 2; ~ 2, while most of the sources with S350/S250 > 0.85 
lie at 2 < 2 < 3 (Amblard et al. 2010; Cooray et al. 2010). 

Finally a bright star mask is applied to all samples 
described above. We follow the procedures in Waddington 
et al. (2007) and mask a circle around all K ^ 12 point 
sources in the 2MASS catalogue within a radius R given by 
log 7i(arcsec)=3.1 — 0.16K. This radius is more conserva- 
tive compared to the star mask used in the public release 
of SWIRE catalogues. In Table [T] we list the number of 
sources, the median redshift and the selection criteria for 
the foreground and background samples. Fig. [1] shows the 
N{z) for each sample. The N{z) of the background is our 
biggest source of uncertainty. If it is a good approximation, 
then B2 is almost completely separated from the foreground 



samples, while Bl has a small overlap with the foreground, 
in which case Wcc is non-zero. 



4 MEASURING THE CROSS-CORRELATION 
SIGNAL 

The cross-correlation between populations 1 and 2 is the 
fractional excess in the probability relative to a random dis- 
tribution (Peebles 1980). We use a modified version of the 
Landy-Szalay estimator (Landy & Szalay 1993) to measure 
the angular cross-correlation function, 

D1D2-D1R2-D2R1+R1R2 

■uncross (e) = , (12) 

where D1D2, D1R2, D2R1 and R1R2 are the normalised 
datal-data2, datal-random2, data2-randoml and randoml- 
random2 pair counts in a given separation bin (see Blake et 
al. 2005 for a disccussion of different estimators of WcroBs)- 
For the foreground samples, we generate random catalogues 
by distributing sources using a uniform distribution. It is 
more complicated to generate random catalogues for the 
background samples. To take into account the noise prop- 
erties in the sub-mm maps and the angular resolution of 
SPIRE, we make maps of randomly distributed sources 
which are processed by the SPIRE Photometer Simulator 
(SPS; Sibthorpe et al. 2009) for observational programmes 
exactly the same as the real data. The catalogues ex- 
tracted from the SPS simulations are then used as random 
catalogues. To reduce shot-noise in the data-random and 
random-random pair counts, our random catalogues (after 
applying the bright star mask) contain roughly 10 times 
more sources than the real catalogues. We use 40 bootstrap 
realisations of the foreground and background samples to 
estimate the errors and covariance matrix. 

As described in Section 2, we need the bias factors 
of the foreground and background samples to calculate 
the expected clustering-induced and lensing-induced cross- 
correlations. In the past, sub-mm sources have been shown 
to cluster strongly (Scott et al. 2002, 2006; Blain et al. 2004; 
Farrah et al. 2006; Blake et al. 2006; Viero et al. 2009). 
More recently, the linear bias factor has been measured to 
be 3.2 ± 0.5 for sources with S'350 > 30 mJy and 3.4 ± 0.6 
for sources with S350/S250 > 0.85 (Cooray et al. 2010). To 
derive the bias factors of the foreground samples, we esti- 
mate the angular auto-correlation function of Fl and F2, 
which can be described by a power-law Wauto = AO The 
amplitude of lUauto is related to the correlation length of the 
spatial correlation function ^(r) — (r/ro)"'''"'""'^' (e.g. Efs- 
tathiou 1991), 



A = frZ / x'-'' {N{z)YE{z)dz[ / N{z)dz 



(13) 



where / = V^V[{'y - l)/2]/r(7/2), E{z) = {Ho / c)[^r,-,{l + 
zf + Q.aY''^ and we have assumed constant clustering in 
comoving units. Finally, we derive the linear bias factor of 
the foreground using the dark matter correlation function 
6 = K(ro)/Cdm(ro)]^/^. The linear bias factor of Fl and F2 
derived in this way is ~1.5 and 1.6 respectively. 

The measured angular cross-correlations between the 
various foreground and background samples are shown in 
Fig. [2I A set of logarithmically spaced angular separation 
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bins are used, ranging from ~ 1 to 50 arcmins. The green 
dashed hne is the expected lensing-induced cross-correlation 
Wfb(S), the red dashed line is the expected clustering- 
induced cross-correlation Wcc{9) and the blue dashed line 
is the sum of the two. In the left column of Fig. [2l the ex- 
pected clustering-induced cross-correlation lOcc is non-zero 
because the tail of the background N{z) overlaps slightly 
with that of the foreground N{z). Although Wcc is much 
smaller than wa,, we should bear in mind that Wcc could be 
underestimated if a higher than expected fraction of SMGs 
reside at low redshifts z < 1. In the right column of Fig. [2l 
the predicted Wcc vanishes, as B2 does not overlap with Fl 
or F2. To assess the significance of the lensing-induced cross- 
correlation signal, given the covariance matrix obtained from 
bootstrap realisations, we derive the Bayes factor 



K = 



P(J|Mlensing) 
P(D|Mn,ll) 



(14) 



where P(-D|Micnsing) is the probability of the data given 
the lensing model and P(_D|Mnuii) is the probability of the 
data assuming there is no cross-correlation. We find that 
K = 6.3 for the cross-correlation between Fl and B2 and 
K = 132.6 between F2 and B2. On Jeffreys' scale (Jeffreys 
1961), K > 3 means that there is substantial evidence that 
A^icnsing is more strongly supported by the data than the 
null hypothesis and K > 100 means that there is decisive 
evidence that Miensing is the favoured model compared to 
the null. Note that there is almost a factor of two increase 
in the source density in the foreground sample F2 compared 
to Fl; increasing the number of tracers of the foreground 
structure increases the strength of the lensing signal. 



5 THE EFFECT OF WEAK LENSING ON THE 
NUMBER COUNT OF SUB-MM SOURCES 

The effect of lensing on the number count of the sub-mm 
sources is expressed in equation (2), under the assumption 
that the lens plane is at a much lower redshift than the 
source plane. The power-law slope of the intrinsic / unlensed 
number count Nu{S) is not affected because the lensing mag- 
nification fl is independent of the flux density. However, the 
overall normalisation of the number count can be modified 
by a factor of /i''"^, where /i = 1 + Sfi — 1 + 2k, in the weak 
lensing limit. Weak lensing by large-scale structure causes 
Sjj, to follow a Gaussian function with mean magnification 
(Sfi) = and its dispersion tj^ dependent on the redshift of 
the sub-mm population (BSOl). Therefore, when averaged 
over a statistically representative area, the effect of weak 
lensing on the number count should be negligible. 

The effect of weak lensing on the local number density 
of the sub-mm sources along a certain direction can be esti- 
mated from the measured cross-correlation between the fore- 
ground and the background populations. In the right panel 
of Fig. 2 where the measured signal is expected to be due 
to lensing only, we can see that the probability of finding a 
background sub-mm source close to a foreground galaxy is 
increased by a few percent above random on angular scales 
between ~ 1 and 50 arcmin. Therefore, the lensing induced 
change in the number density along a certain direction is 
expected to be at the level of a few percent. 

We can also estimate the effect of lensing on the local 



number density through the auto-correlation function of the 
background sub-mm sources, TOauto(^) = {&n{(f))5n{(l)')) . Us- 
ing equation (4), we can decompose uiauto(^) into three com- 
ponents, {5rf{^)Sn''(4)')), {5Tf{^)5n^'{4)')) + {Sn''{<j)')Sn''{ij))) 
and {5n'^{(j>)5n'^ {(!>')), which represent the galaxy- galaxy, 
galaxy-lensing and lensing-lensing correlation functions re- 
spectively. The lensing-lensing term is given by (Moessner 
& Jain 1998) 



lensing- 



-lc„sing(g) = [3an(/3-l)]' / {g2/afdx 

Jo 



^P(x,k)Mkr9)dk. (15) 

At zero lag, wiOyZT"-''""""" = {(Sn^f) is the variance of 
the number density fluctuatuation due to lensing and thus 
the rms fluctuation is 5n^' = («^(0)Luto"''"'""""'')''" which is 
at a few percent level. 



6 DISCUSSIONS AND CONCLUSIONS 

The unusually steep number count in the bright sub- 
mm regime leads to an enhanced cross-correlation signal 
that is due to weak gravitational lensing. In this paper, 
we have measured the angular cross-correlations between 
sub-mm sources detected by _ffersc/ie/-SPIRE in Lockman- 
SWIRE and foreground sources selected in the optical or 
near-infrared. We have also derived theoretical expecta- 
tions of the weak lensing-induced cross-correlation liifb and 
the clustering-induced cross-correlation Wcc which are in 
good agreement with our measurements. We find clear ev- 
idence for a lensing-induced cross-correlation between sub- 
mm sources at high redshifts and galaxies at low redshifts. 

The redshift distribution of the sub-mm sources is the 
biggest source of uncertainty in our analysis because most of 
the sources do not have spectroscopic redshifts. In principle, 
the clustering-induced cross-correlation Wcc could contam- 
inate the lensing-induced cross-correlation uifb if a higher 
than expected fraction of sub-mm sources reside in the low- 
redshift Universe. As the amplitude of uifb is mainly sensitive 
to the mean redshift of the background population rather 
than the exact shape of the N{z) (Menard & Bartelmann 
2002), we have carried out a simple calculation of the ex- 
pected uifb and uicc amplitude by varying the mean redshift 
{z) (from 0.3 to 4.0) and the width (from 0.2 to 2.5), 
assuming the N{z) of the sub-mm sources can be approx- 
imated by a Gaussian function. In all cases, to reproduce 
the measured cross-correlation signal, Wcc is at most com- 
parable to uifb when (z) ~ 3.5, ~ 1.5, (z) ~ 2.5, ~ 1.0 
or (z) ~ 1.5, ~ 0.5. So the detection of the weak lensing- 
induced cross-correlation should be robust. It should be pos- 
sible to acurately determine N{z) in the future when the 
infrared spectral energy distributions are well understood 
and/or more spectroscopic redshifts are acquired for sub- 
mm sources. 

Limitations in our modelling of the cross-correlation in- 
clude: using a scale- and time-independent bias factor for the 
galaxy-dark matter power spectrum; assuming a linearised 
magnification; and adopting a constant power-law number 
count slope independent of flux. While for this first study 
a simple model is adequate given the large error bars, an 
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Figure 2. The angular cross-correlations between foreground and background populations. The error bars are the rms scatter derived 
from 40 bootstrap realisations of the real data. In each panel, the red dashed line is the predicted cross-correlation due to mutual 
clustering Wcc{6)- The green dashed line is the predicted cross-correlation due to lensing W{^,{6). The blue dashed lines show the sum of 
Wcc{S) and W{^,{6). The black horizontal line denotes the zero level. In the right column, the expected Wcc{9) = 0. 



approach such as the halo model to describe the galaxy- 
dark matter power spectrum can be utilised in the future 
when additional data warrant an improved description (e.g. 
Jain et al. 2003). The expected increase in area covered by 
Herschel-SPIRE will allow the detection of cosmic magnifi- 
cation presented in this paper to be improved and be used 
to constrain cosmological parameters and galaxy bias. 
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